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Abstract 

We discuss hadronic effects in the muonium hyperfine structure and derive an 
expression for the hadronic contribution to the hfs interval in form of the one- 
dimensional integral of the cross section of e + e~ annihilation into hadrons. Higher- 
order hadronic contributions are also considered. 
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1 Introduction and Results 



Quantum Electrodynamics (QED) provides an opportunity to determine char- 
acteristics of various particles and simple atomic systems. However, any pure 
QED caclulation is incomplete even in the case of a purely leptonic system 
because of a contribution of the strong interaction originating from hadronic 
intermediate states. For example, a contribution of the hadronic effects to the 
anomalous magnetic moment of the muon (a M ) is about 60 ppm of its QED 
value. This value is larger than both the uncertainty of the QED calculations 
(see e. g. [1]) and experiment [2]. The leading part of this correction (Fig. 1) 
can be presented in the form 

Acv(hadr) = ^- 2 j^-K a {s) R{s) . (1) 

Here a is the fine structure constant and relativistic units in which h — c—1 
will be used throughout the paper. Such a presentation is quite useful since it 
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Fig. 1. The hadronic con- Fig. 2. The hadronic contribution to the muo- 

tribution to the anomalous nium hyperfine splitting 
magnetic moment of muon 

clearly separates the QED and hadronic parts. The QED part is known [3] in 
a closed analytic form 
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while the hadronic factor 

_ a(e + e- -» 7 -> hadrons) 



can be determined either from direct measurements of the cross sections of 
e + e~-annihilation into hadrons at the energy E cra ^ = ^/s or from theoretical 
estimations. Here r = (1 — Am 2 J s) ' . 

In this paper we derive a similar expression for the hyperfine structure (hfs) 
interval in a two-body hydrogen-like atom with a point-like nucleus (Fig. 2). 
In particular, in the case of the ground state hfs in muonium the result is 
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where 



^-W - - (£j +2) r h i±I + (^ + I) in^ - I , (4) 



and £jr = 8a 4 m|,/3m e m M is the leading contribution to the hyperfine struc- 
ture (the so called Fermi energy) and ttir is the reduced mass. 

Presently, the accuracy of the QED calculations [4-6] and of evaluations of 
the hadronic vacuum polarization contributions (see e.g. [7]) is sufficient for 
comparison with the experiment [8]. However, the expected progress in de- 
velopment of intensive muon sources for needs of particle physics offers an 



opportunity to increase the statistics of muonic events and to provide a much 
better source of muonium [9]. In view of this increase of the experimental ac- 
curacy we need more precise theoretical calculations of various contributions 
and in particular of those for the hadronic vacuum polarization (Fig. 2). The 
accuracy of such a hadronic calculation should establish the frontiers of any 
possible precision test of the bound state QED with muonium. 

The expression (3) was obtained by us some time ago and the result appeared 
in a paper of one of us [10]. However, neither derivation nor discussion of the 
corrections to (3) was presented. Since a calculation based on this expression 
is now in progress [6], the derivation of (3) will be given here in detail. The 
accuracy of an incoming calculation [6] based on our expression (3) is of the 
1%-level and we discuss here various corrections to (3) and in particular those 
due to the higher-order hadronic vacuum effects (Fig. 3). Eventually we found 
that there is only one higher-order correction to (3) above the uncertainty level 
of the calculation in [6] and it increases the vacuum polarization contribution 
by roughly 3%. 
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Fig. 3. Second-order hadronic corrections to the muonium hfs 



2 General expression 



We start the derivation with the two-photon exchange diagrams. Their con- 
tribution can be presented in the form (cf. [11]) 



= g ^ 



>Sc 



7T m; - m% 



I{m^) - I(m e 



where Km) 



dk 2 



J(k,m) 



and 



J(k,m) 



1 



2 ( -Vk 2 + Am 2 - 1 
,k 



Am 2 



kVk 2 + 4m 2 -k 2 -2m 2 ) . (5) 



When calculating the integral over the euclidean momentum k a divergency at 
low momentum should appear and it is necessary to rearrange the expressions. 
By introducing the vacuum polarization insertion in the skeleton integral (5), 
we reach our goal without any additional rearrangement. 

We take the polarization into account substituting the photon propagator with 
a dispersion presentation of the polarization 
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(6) 



The choice of the dispersion weight function p(s) and the limits of the inte- 
gration depend on what contribution is to be calculated. In the case of the 
hadronic vacuum polarization one finds s > 4m 2 ^> ml, and p(s) = R(s)/3s 
where R(s) is defined in (2). Finally, we present the vacuum polarization con- 
tribution in the form 
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(7) 



where 



I VP {s,m) = I 



dk 2 



k 2 + s 



J(k,m) 



After integrating over the momentum k we find 
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That is the most general expression for the hadronic contribution to the hy- 
perfine splitting in a hydrogen-like system with a point-like nucleus. In the 
case of the hadronic contribution to the hfs interval in muonium one can take 
advantage of s > ^m 2 3> 4m 2 and neglect a contribution of m — m e 
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Fig. 4. The QED kernels for the hadronic contributions to the muonium hyperfine 
splitting (Kyiu) and to the anomalous magnetic moment of muon (K a ) and their 
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leading asymptotics, K^ lu = (181n(s/ra„) + 15)m^/4s and K a = m^/3s. The 
kernels are shown in Fig. a with solid lines, while the asymptotics are with dashed 
lines. In Fig. b a result for (Kmu — ^Mu)/^Mu is presented. 

After simple transformations we arrive to the result (3). This expression is 
appropriate for calculations of any vacuum polarization contributions with 
s < 4:m 2 such as the hadronic or r-leptonic contributions. The neglected term 
(9) does not exceed 10" 3 of the main contribution for s > (2m 7T ) 2 . 



3 Asymptotics of the QED kernels Kmu(s) at high energy 



Since s > (2m 7T ) 2 > (2m M ) 2 the asymptotic behaviours of the QED kernel 
K Mu at high values of s is of interest. We find 
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(10) 



In the case of the p-meson contribution, which is dominant (see Sect. 4), the 
parameter of expansion is 4m 2 /m 2 ~ 0.076 and the leading contribution at 
s = m 2 is 0.409, while the higher order terms are only 0.006. 

This behaviour should be compared with the expansion of K a (s) (see Fig. 4) 
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4 Estimation within a simplified narrow-pole model 



The essential part of the hadronic vacuum polarization contribution comes 
from the p- meson-pole and it can be easily estimated (cf. [12]) from a simple 
approximation 
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with a sum over three light vector mesons dominating in e + e~ annihilation 
at low energies: p, u and 0. The partial width of a resonance decay into 
an electron-positron pair is a value determined experimentally. The results for 
the p-, io- and 0-mesons within the simple pole approximation are summarized 
in Table 1 as well as their masses and electronic widths. This calculation is 
important to check if our expression is consistent with previous calculations of 
the hadronic contributions. The p-meson contribution is 65% of the hadronic 
contribution [6] and all three pole contributions deliver 76% of Az/hadr (or 
0.18 kHz or 4.0 • 10~ 8 z/p). The value of the p-meson contribution is in a good 
agreement with the 7T7T contribution in [6] . Due to that we study any correction 
to (3) with the help of these resonance contributions obtained from the simple 
pole approximation. 



Resonance 


p-meson 


w-meson 


0-meson 


M res [MeV] 


769(1) 


783 


1019 


T(res — ► e + e~) [keV] 


6.77(32) 


0.60(2) 


1.32(5) 


Az/ rcs [kHz] 


0.151(7) 


0.013 


0.014 


Afrcs/Az^hadr 


65(3)% 


5.5(2)% 


6.0(2)% 


-^Mu 


0.415 


0.404 


0.262 


A Mu 


0.409 


0.398 


0.260 



Table 1 

Vector meson resonances, their parameters and contributions to the hyperfine struc- 
ture. The fractional value in respect the to the complete hadronic contribution is 
given for Ai^hadr = 0.233(3) kHz [6]. The parameters are taken from [13] 



5 Higher-Order Hadronic Contributions 



To estimate the higher-order hadronic contribution in order a 2 (Za) we simply 
note that only one set of diagrams leads to a contribition enhanced by a 
logarithmic factor Yn.(s/m 2 e ) (Fig. 3a). These diagrams contain the electronic 



vacuum polarization insertion and in the logarithmic approximation can be 
easily calculated for the p-meson contribution: 
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The same estimation can be reproduced for two other resonances in Table 1. 
About 80% of the hadronic contribution comes from these three resonances, 
while the rest is essentially related to the high-s backdround. It it easy to check 
that in the case of the background contribution and other resonances the loga- 
rithmic approximation is valid and we confirm the logarithmic approximation 
(cf. (13)) for the complete result, which is about 3% of the leading contri- 
bution. We expect that a non-logarithmic part of the cotributions (Fig. 3) is 
below 1%-level. 

A similar estimation is misleading in the case of the anomalous magnetic 
moment of the muon [14], where the contribution related to insertion of the 
electronic vacuum polarization into leading hadronic diagram Fig. 1 does not 
dominate among higher-order hadronic contributions. We believe that it is 
caused by a special structure of the kernel (11), which has a small numerical 
coefficient and no logarithmic enhancement at high s. That is not a case 
for the muonium hfs (10) and we finally estimate the higher order hadronic 
contribution as (3 ± 1)% of the leading hadronic term, or 0.007(2) kHz. 
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